Generation of entanglement density within a reservoir 
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We study a single two-level atom interacting with a reservoir of modes defined by its reservoir 
structure function. Within this framework we are able to define a density of entanglement involving 
a continuum of reservoir modes. The density of entanglement is derived for a system with a single 
excitation by taking a limit of the global entanglement. Utilizing the density of entanglement we 
quantify the entanglement between the atom and the modes, and also between the reservoir modes 
themselves. 
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I. INTRODUCTION 

In recent years, entanglement has attracted the atten- 
tion of many physicists working in the area of quantum 
mechanics [l|, |2| . This is due to the ongoing research in 
the area of quantum information Q , and also because of 
the advances made in different experimental disciplines, 
such as in ion traps Q and Bose-Einstein condensation 
0, Q. Developments in the field of cavity QED, where 
experiments in the strong coupling regime are carried out 
[2, iS] provide plenty of motivation for studying quantum 
information and entanglement. Theoretical studies are 
also important in the context of atom-light interactions 
inside structured reservoirs Q such as resonant cavities 
or photonic band gap materials. The theoretically pre- 
dicted atom-photon bound state could also lead to entan- 
glement and this can also be linked to another problem: 
that of atom-laser out-coupling from Bose-Einstein con- 
densates (lol - [r3 |. where analogous effects were predicted 
in the past. 

In a recent work by Cummings and Hu 13], a two- 
level atom coupled to a large reservoir was considered. 
Starting with simple models with a few modes, they gen- 
eralize to reservoirs with a large number of modes, where 
they examine entanglement between the atom and the 
reservoir. In their analysis the reservoir is treated as a 
collective object, i.e. the entanglement between the atom 
and each mode, or between the modes, is not considered. 

In the present work we study entanglement between 
the atom and the reservoir by means of a different ap- 
proach. In our model, an atom interacts with a contin- 
uum of modes at zero temperature [3, 113 , arid entan- 
glement properties between the atom and the modes and 
also between the modes are considered. Using global en- 
tanglement as a measure of entanglement, we derive a 
pair of distributions that can be interpreted as densities 
of entanglement in terms of all the reservoir modes. Both 
distributions can be calculated in terms of the spectrum 
of reservoir excitation. With these two new measures of 
entanglement we can study in detail entanglement be- 
tween the atom and the modes, and also between the 



modes. 

In our analysis, we consider a Lorentzian reservoir and 
cover different dynamic regimes. For strong coupling we 
observe the periodic collapses and revivals of entangle- 
ment between the atom and the reservoir, which are asso- 
ciated with Rabi oscillations of the atom. Eventually, all 
the population leaves the atom and the reservoir becomes 
entangled. More precisely, two bundles of modes are ex- 
cited, forming a "Bell-like" state. The method developed 
here, in terms of entanglement distributions, can also be 
used when considering various types of structured reser- 
voirs. For example, reservoirs with a density of modes 
characterized by a band gap d, [H, [l3| can be treated 
with this method. 

The paper is organized as follows. In section |TT] we 
introduce a measure of entanglement, the global entan- 
glement, and apply it to our system-reservoir states when 
the number of degrees of freedom is finite. In section Hill 
we introduce the density of entanglement for the limit of 
a bath with an infinite number of modes. This entangle- 
ment measure evolves according to our model system as 
introduced in section llVl In section|V]we present our key 
results for the density of entanglement and we conclude 
in section IVll 



II. MEASURES OF MULTI-PARTITE 
ENTANGLEMENT 



Identifying and measuring entanglement in multi- 
partite systems presents various complications. Apart 
from the case of a two-qubit system, where entangle- 
ment can be identified both for a pure and a mixed state 
[l9| . multi-qubit entanglement is an open problem 
and to date several measures of entanglement have been 
proposed P, 0, [20l - [23l | . For the analysis that follows we 
will be using the global entanglement [l^l since this will 
enable us to deal with many modes (or qubits). The 
physical problem we consider is that of entanglement be- 
tween an atom and a reservoir, and in the context of 
this problem we will also consider the entanglement be- 
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tween different reservoir modes. To quantify tfiis we will 
take a discrete bath of reservoir modes. Since we will 
assume only a single excitation in each bath mode, we 
can treat the bath states as a set of qubits for the pur- 
pose of computing the entanglement, and then later take 
a continuum limit. With just one excitation in total in 
the system, this excitation may be in the atom, or in the 
bath, so that the state vector at all times will have the 
form 

IV^W) =Co|0) +Ca(i)|V'a) +5ZcaWIV'a), (1) 

A 

where |0) is the vacuum state 

|0) = |0)a(8) 10000 •••000). (2) 

The other two states |V'a) and |V'a), correspond to the 
atom being excited 

|^a) = |l)a®|000^^^000), (3) 

or the mode A of the reservoir being excited 

I^a) = |0)a® IOOO^^^OIaO^^^OOO). (4) 

Recently this kind of approach has been utilized for the 
study of decoherence and entanglement decay in systems 
with one or two excitations (effectively at T = 0) [23 - 
1131 . Some of this work examines the decay of a single 
excitation in a reservoir made from a continuum of modes 
[l^ll^jllOjIsO] and other works examine entaglement with 
two or more qubits, but most still utilize the decay of a 
single excitation as considered in the present paper [23 - 

[13, iSii EHil Ei, S, Ei]. 

If we started with a system of N qubits in a pure state 
IV'), the global entanglement is defined as 

9 ^ 

Q(|V'))=2--^trp?, (5) 

where pi is the reduced density matrix for the i-th qubit. 
When is a pure product state then Q{\ip)) = 0. If we 
take an entangled state for such as the GHZ state 

\GHZ)n = (|000 • • • ) + |111 • • • ))/%/2, (6) 

then Qdip)) — 1. This example gives a maximum value 
for the global entanglement which is normalized such that 
< Qdij})) < 1. Another, more relevant example is the 
W-state 

N 

\W)n = Y.\0---01,0---0}/Vn, (7) 

for which the global entanglement goes to zero as 1/N 
for large N 

Q(|V)) =4(iV-l)/iV^ (8) 



An important property of the global entanglement is 
that it is equal to a sum over two-qubit concurrences . 
More specifically, for pure states \tp{t)) [Eq. ([1])] we find 
from Eq. ([5]), see the appendix for details, that 

QM) = (9) 

where C^(i) reads 

N 

A=l l<\<fi<N 

The concurrence c^{pa\) is that for the two-qubit (re- 
duced) density matrix pa\ [H, [111 

Pax=tr,^)^{\m)m)\) (11) 

between the atom and the A-mode. The quantity p^A is 
the corresponding density matrix for the modes pL and A: 

P^x = ira,.^x..A\4'{t)){m\}- (12) 
The two-qubit concurrence c{p) is equal to [H, [l^ 

c{p) — max{0, — \/\2 — ~ V^}, (13) 
where \j are the eigenvalues of the matrix 

R = p{(Ty® <Jy)p*{ay® Uy), (14) 

in decreasing order, i.e. Ai > A2 > A3 > A4, and Uy is 
the relevant Pauli matrix. 

Is is important to note here that (t) is exactly equal 
to the square of the norm of the concurrence vector [23| , 
which is one of the many proposed measures for multi- 
partite entanglement. Because of the connection between 
C^{t) and the two-qubit concurrence Eq. (fTO|l . we shall 
refer to C^{t) simply as the concurrence. 

For the remainder of this work we will focus only on 
C^{t) and its properties since, with the exception of the 
normalization factor 2/{N -I- 1) in Eq. ([9]), it is equiva- 
lent to the global entanglement. Furthermore, we will be 
considering a continuum of reservoir modes, i.e. A'^ ^ 1. 
In this limit, as we see in the following section, one can 
define a density of entanglement for continuous systems. 

III. DENSITY OF ENTANGLEMENT 

In the limit of continuous distribution for the reservoir 
modes, the summations in Eq. (fTO|) can be converted into 
integrals over the density of modes p\, i.e. 

E ^ / dc^APA. (15) 

A 

If we take this limit, the concurrence C^(t) can be ex- 
pressed as the sum of two separate parts which ultimately 
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involve either atomic population or reservoir mode pop- 
ulations: 



(16) 



OO pOO 



+ I dojx / dLLii_,£R{uj\,ujfj.,t). 

-OO J —OO 



The two distributions £a and En will form the entangle- 
ment densities and are defined in terms of the two qubit 
concurrences to be 



£A{iO\,t) ^ C^{pa\)p{l^\)., (17) 



and 



£RiuJx,UJf„t) = -C^(pXf,)p{uJx)p{ujf,), (18) 

where p{u!) is the reservoir density of modes. The inter- 
pretation for these two functions is now very simple. 

The distribution SAii^Xjt) is the density of entan- 
glement between the atom and all reservoir modes in 
the vicinity of mode uj\ i.e. £yi('^A, i)dwA is the to- 
tal entanglement between the atom and all modes in 
the frequency interval ujx and uj\ + duj\. In the same 
way £ji{uj\,ujfj_,t)duj^ gives the entanglement between 
the mode wa and modes in the frequency range uj^ to 
-I- dw^. 

For the remainder of this work, £a will be referred 
as the atom-mode density of entanglement and to £fj 
as the mode-mode density of entanglement. The global 
entanglement Qdtp)), Eq. ([9]), can be calculated from 
these two densities of entanglement, and thus a number of 
entanglement properties can be studied in terms of these 
two distributions. To this end it is important to note an 
interesting property for the two distributions £a and Sr. 
One can show that both entanglement distributions can 
be written in terms of the spectrum of reservoir excitation 

S{LO,,t)=p{LOx)\cx{t)\', (19) 



as 



and 



£AiuJx,t)^A\cait)\^S{uJx,t), (20) 



£R{ux,uj^,t) = 2S{ux,t)S{iu^,t) . (21) 

Thus, both entanglement distributions can be derived 
from the reservoir excitation spectrum. Having the def- 
inition for the entanglement density we can now intro- 
duce a Hamiltonian and dynamics to consider the time- 
dependent properties of entanglement for an atom cou- 
pled to a Lorentzian reservoir. 



IV. DYNAMICAL MODEL 

The model system consists of a two-level atom coupled 
to a reservoir of harmonic oscillators with annihilation 
and creation operators a\ and respectively. Within 
the rotating wave approximation the Hamiltonian reads 



H =^^uJxa\ax +a;o|l)aa(l| 
A 

+ E-9^ (4|0)„„(l|+aA|i; 



(22) 



where gx is the coupling between the mode A and the 
atomic transition |l)a — > \Q)a- The atomic transition 
frequency is ujq whereas the A-mode frequency is uJx- This 
model will preserve the assumption of a single excitation 
which is built into the system states in Eq. 

For the purposes of the analysis that follows, it is very 
useful to introduce the reservoir structure function Dioox) 
which reflects the properties of the density of modes [ig] . 
This is defined through 



px{gx? = ^D{u.x\ 



and is normalized such that 



/OO 
dujD{Lj) ^ 2tt. 
-OO 



(23) 



(24) 



With this normalization a measure of the overall coupling 
strength is flo which is given by 



(25) 



Utilizing these assumptions, and a state vector of the 
form given by Eq. ([T]), the Schrodinger equation in an 
interaction picture yields 



(26) 



(27) 



with the detuning between the atomic transition and the 
mode A being given by 



6x=ujx- ujo- 



(28) 



The interaction picture amplitudes Ca and cx are con- 
nected to Ca and cx via a time-dependent transformation 

Ca{t) = e'""*c,(t), dxit) = e'"-*CA(t). (29) 

Apart from numerical integration of equations (j26p and 
(P7)l . one could use analytical methods to derive the dy- 
namics. Examples of these are: the resolvent method [9j, 
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the Laplace transform [56|, |57[ and that of the pseudo- 
modes This latter method applies when 
the spectral function is analytic with poles in the lower 
complex plane. Then equations (pS)) and (P7|) can be 
replaced, in the continuum limit, by a set of equivalent 
equations where the atom now couples to a finite set of 
fictitious modes. Each of these modes has a one-to-one 
correspondence to the poles of D{uj). 

For concreteness we consider Lorentzian reservoir 
structure functions of width F and peak frequency Wci 
i.e. 

Diu;) ^ 



(L.-c.,)2 + (r/2)2- 

This reservoir function has a simple pole at 

z ^ ujc- iF/2, 

and for this D{uj) equations p6| and 
to two new equations 

. d , , , 



dt 



(A - iF/2)6(t) 



where the atom-cavity detuning A 



bit) 



-iujQt 



h{t) is the pseudomode amplitude. 



(30) 



(31) 

27|) are equivalent 



(32) 

(33) 

Wo , and 
This set 



of equations can be associated with a master equation, 
where the physical interpretation for the pseudomode is 
that of a leaking cavity mode coupled to the atomic tran- 
sition 16]. 

These coupled ODEs are straightforward to solve, and 
in particular, for the case of resonance, Wc = wq, one finds 
for the atomic amplitude Ca(t): 



Ca{t) = Ca{Q)e'' 



cosh {at 1 4) ^ sinh (at/4) 

a 



(34) 



The modified decay rate a is given by 



(35) 

in terms of F and VLq. For strong couplings, 4r2o/r ^ 1, 
the atomic population oscillates between the atom and 
the reservoir with a slowly decaying amplitude. In the 
long time limit all the energy is lost to the reservoir. 
If we decrease the coupling rig, the atom dissipates its 
energy faster, and for 4fio/F ^ 1 the atom exponentially 
decays into the reservoir. These two types of behavior are 
illustrated in Fig. [1] which is given as a reference point 
for the entanglement discussion in the next section. The 
figures show the atomic population for 51o — lOF and 
ilo = O.IF plotted against time t. 



V. EVOLUTION OF RESERVOIR 
ENTANGLEMENT 

A. Entanglement generation by decay 

The interaction between the atom and the reservoir re- 
sults in entanglement creation between the atom and the 




FIG. 1. The atomic population |ca(t)|'^ as a function of time 
for strong coupling f2o = lOF (solid line), and weak coupling 
Slo = O.ir (dashed line). Both results were obtained from 
Eq. (|34p . The initial system state in this and the following 
figures is defined by Ca(0) = 1, cq = and all Ca(0) = 0. The 
detuning A = 0. 



reservoir modes, and also between the reservoir modes. 
This latter entanglement is indirect and is due to the ef- 
fective coupling between the modes as a result of their 
interaction with the atom. When considering our mea- 
sure of the total entanglement, the concurrence C^it) as 
a function of time (Fig. [5]), we see that for both strong 
and weak coupling the concurrence builds up to a maxi- 
mum value equal to C^(i) = 2. In the case of strong cou- 
pling, Fig. 2(a) the concurrence reaches the maximum 



value very quickly, and then oscillates in a way which 
is closely connected to the Rabi oscillations seen in Fig. 
[H The oscillations return to the maximum value of two, 
with a minimum value that also approaches two as time 
increases. If we compare the concurrence to the popula- 
tion in Fig. [T]we see that every minimum in concurrence 
is matched by a peak in population of the atomic state. 
Likewise the maxima in concurrence are matched by min- 
ima in the atomic population: it seems that the decaying 
atom is very efficient at generating entanglement in the 
reservoir. At this point we note that at long times the en- 
ergy has left the atom (see the population Fig.[T]) and the 
system is in an approximate product state of unexcited 
atom and bath states. The large value of the concurrence 
at these times indicates the presence of entanglement in 
the bath. In the long time limit the concurrence reaches 
a steady value C^(oo) = 2, a result that can be derived 
analytically if one uses Eqs. (fT6)) pO)) and (l2Tt and the 
spectrum for f — > oo. This limit for the concurrence is 
the same in the weak coupling case, see Fig. |2(b)[ the 
atom decays exponentially and the entanglement reaches 
a steady state monotonically. In the long time limit the 
atom again disentangles from the reservoir and entangle- 
ment is distributed only between the reservoir modes. 
By examining the entanglement densities [Eqs. (|17I18I) ] 
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(a) 




5 10 15 20 



£lot 

(b) 

FIG. 2. The concurrence C^{t) as a function of time for dif- 
ferent couplings 0,0 = lOr (a) and Qo = O.IF (b). The resuh 
was obtained from Eqs. (fT6)) . JSOj and pT]) . 




FIG. 3. The density of entanglement f a(wa, i) between atom 
and bath is plotted for different times t and reservoir mode 
frequencies The coupling strength is: (a) fio = lOF; and 
(b) flo = O.IF. 



we can gain further insight into where the entanglement 
resides in this system and how it evolves over time. In 
Fig. [3] we show the atom-bath density of entanglement 
£a{uj\, t) for both strong and weak couplings. The strong 
coupling case shows a complex behavior. First a central 
peak of entanglement appears in the vicinity of uj\ ~ wq; 
note the peak at (5a ~ for short times in Fig. 3(a) The 



entanglement is then transferred to Rabi sidebands at 
d\ = 5± ^ zL^Iq. The sideband entanglement oscillations 
seen take place at half the frequency of the central peak 
entanglement oscillations. Ultimately all the entangle- 
ment decays at long times, £A{tox,t) — )■ 0, as there can 
be no entanglement between bath and atom when the 
atomic population approaches zero. At that point the 



entanglement indicated by the concurrence in Fig. 2(a) 



must reside in reservoir mode entanglement which is ex- 
amined in the next section, Sec. lVBl Fig. 3(a) shows that 
as Rabi sidebands develop in the reservoir excitation, the 



atom-bath entanglement moves from central frequencies 
to the sidebands at 5\ = 6±. The period doubling of the 
central peak is due to the oscillations of the excitation 
there, combined with the oscillations of atomic popu- 
lation. The population of the sideband modes is more 
stable, but oscillations of the atomic population result in 
oscillations of entanglement there, too. 



For weak coupling, Fig. 3(b) the entanglement briefly 
resides across the whole reservoir structure. This is es- 
sentially because the coupling of the atom is over this 
same range. In weak coupling, however, Rabi sidebands 
do not develop. Instead the final population of bath 
modes will be over a relatively narrow frequency range 
well within the model Lorentzian profile, Eq. pOI) . Since 
both atomic, and mode population is needed for atom- 
bath entanglement, the narrow central frequency region 
is entangled only for a short while and then decays. 
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B. Entanglement between reservoir modes 



As already mentioned, coupling the atom and the 
reservoir modes will induce an indirect coupling between 
the reservoir modes. Because of this, the modes will en- 
tangle and it is important to consider the properties of 
the mode- mode density of entanglement Er. In the long 
time limit one can easily obtain analytic expressions for 
the spectrum of excitation in the reservoir and from that 
calculate the density of entanglement for i — >■ oo. Using 
the definition for the spectrum of reservoir excitation Eq. 
(fTO]) with the solution ^ and Eq. we find that for 
t — oo 



2tt 



51 



(r/2)2 



(36) 



VtlT/2 



where, as before, 5x = ujx—ljq [Eq. (|28|) ]. Taking Eq. ([36|) 
together with Eq. (PT|) it is straightforward to calculate 
the mode-mode density of entanglement for i — >■ (X) 

£r{llix,uj^, oo) = 



2n^ 



{T/2y5l 



(37) 



Figure |4] shows the mode-mode density of entangle- 
ment for several couplings and in the limit i — oo. In 
the strong coupling case. Fig. 4(a) we see the formation 



of four sharp peaks which signify the existence of strong 
entanglement between the two symmetric modes 5± as 
a result of the Rabi splitting. This picture has a simple 
and rather intuitive interpretation: the final state for the 
reservoir has the form of a Bell-like state between the two 
symmetric modes 5± w ±fio- First we note that the final 
state takes the approximate form 

N 

IV(cx))) = ® J2 (^+(^^) + ^-(^^)) 1^^)- (38) 

A=l 

where the two probability distributions |P±(ajA)P are 
centered at 6x = 6+ and 6- respectively. For fio ^ T, i.e. 
in the strong coupling regime, their width is very small 
which practically means that only the two modes S± are 
excited and, for this, the reservoir state can be approxi- 
mately described by a Bell state of the form 

)®(|l+0_)+e''^|0+l_))/V2, 

(39) 

where (/) is an arbitrary phase factor. This picture applies 
only at long times. For short times the mode-mode den- 
sity of entanglement will initially have a distribution that 
peaks in the vicinity of ujx ~ uj^ ^ ljq. As time evolves 




(a) 




(b) 




(c) 



FIG. 4. The density of entanglement £R{u)\,ui^,t) between 
bath modes for t — i- oo. The coupling strength is: (a) S7o = 
lOF; (b) fio = T; and (c) fio = O.ir. 
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this distribution breaks into four symmetrical peaks seen 



in Fig. 4(a) as a result of the Rabi splitting. 



no 



For moderate and weaker couplings, ilo = F and 
<C F, the simple patterns observed in Fig. 4(a) 



for 



the mode-mode density of entanglement due to the Rabi 
splitting disappear. For example in Fig. |4(b)] the density 
of entanglement for t — >■ oo and fio = F is plotted. From 
this we can see that, for moderate couplings, although the 



distribution peaks around 



Sx 



6±, the entangle- 



ment spreads over a wider range of reservoir frequencies. 



For even weaker couplings, Fig. 4(c) only modes in the 
vicinity of uj\ « ujq are excited and thus the entanglement 
is created only between modes in this frequency range. 
This is the reason for having a peaked entanglement dis- 



tribution with a centre at uj,. 



be tested with a measure of non-Markovianity |59f . 

In conclusion, when considering entanglement between 
an atom and a large reservoir, the analysis can be formu- 
lated in terms of entanglement density functions. These 
distributions can be associated with the spectrum of 
reservoir excitation, a quantity that can be measured, 
for example, in cavity QED experiments. Although the 
model considered here, an atom coupled to a Lorentzian 
reservoir, is rather simple, it can be extended to consider 
more complicated reservoir structures such as model pho- 
tonic band gaps. Potentially, one could consider general- 
izations for the density of entanglement, and also extend 
the model beyond the assumption of a single excitation 
in the system. 



VI. CONCLUSION 
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In this work we have presented a method for study- 
ing entanglement between a two-level atom and a large 
reservoir of modes such as may be found in, for exam- 
ple, a lossy cavity field. Our aim was to study entan- 
glement between the atom and the reservoir, and within 
the reservoir itself. To do this we restricted our inter- 
est to a simple case where the reservoir modes can be 
treated as qubits. Using global entanglement, we derived 
the entanglement densities £a and S^i for the reservoir 
modes. These distributions have been given in terms 
of the well known two-qubit concurrences, and can be 
calculated from the reservoir excitation spectrum. The 
entanglement densities are then used both for studying 
entanglement between the atom and the reservoir modes, 
and also between the modes. 

In considering different dynamical regimes defined in 
terms of the coupling strength, we noticed that when 
strong interactions occur, the reservoir modes are entan- 
gled in a "Bell-like" state in the long time limit. This long 
time limit is a regime where no excitation remains in the 
atom and all the entanglement is amongst the reservoir 
modes. Since there is no direct interaction between the 
bath-modes, see the Hamiltonian (l22t . the final entangle- 
ment arises through indirect interaction. Another way of 
viewing this is that in the strong coupling regime we have 
a non-Markovian system. In such the atom de- 

cays, some information resides in part of the reservoir in 
a way that it can be returned to the atom later (32j . This 
allows the indirect coupling between the reservoir modes 
which creates the entanglement. In the weak coupling 
(Markovian) case, the information does not return to the 
atomic system and the reservoir entanglement cannot be 
created. 

In the transition from strong to weak coupling, the 
"Bell- like" state of Fig. 4(a) coalesces into a single peaked 
structure as was seen in Fig. 4(c) Based on the interpre- 
tation given above we would expect that the transition to 
a single peak takes place as the system becomes Marko- 
vian rather than non-Markovian. In principle this could 
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Appendix: Global entanglement 

For the completeness of our analysis, we show here how 
the global entanglement for the pure state \tp{t)), Eq. ([T|), 
can be derived. First, the state \ip{t)) can be written as 
a pure state of -t- 1 qubits, i.e. 



m))^co\o)+j2cjit)\t 



3/' 



(A.l) 



where j = a for the atom, or j = A = 1, • • • , A'^ for the 
reservoir modes. The vacuum state |0) is given in Eq. 
([2]), while \tpj) = \ipa) [Eq. (O] for j = a, or is given by 
Eq. (HI) if j = A, i.e. refers to a reservoir mode. 

Next we calculate the reduced two-qubit density ma- 
trix pj, 

Pji = trk^j^iip} = (xoIpIxo) + i^k\p\Xk), (A.2) 
where |xo) is the N — 2-qubit vacuum state 



ixo) = n lOfc) 



(A.3) 



and \xk) is the A^ — 2-qubit state with a single excitation. 



IXfc) = IO---OUO---0). 



(A.4) 
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With these definitions for |xo) and |xfc), the reduced den- 
sity matrix takes the form 



/I 

V 



coc* 
coc* 




I Ci I Cg Ci Cq Cj \ 
12 



CjC* 



C j Ci I Cj I 







(A.5) 



y 



where the basis states, starting from top and moving to 
bottom, are 



|0,0,), |0,1,), |1,0,), 



(A.6) 



The concurrence c{pji), Eq. ()13p . for this density matrix 
reads 



(A.7) 



In order to calculate the global entanglement, we derive 
the single qubit density matrix from Eq. (jA.SP by tracing 
out the ?-qubit 



Pj 



1 - |CjP CqC 



(A., 



which after using probability conservation becomes 



tipj = 1 - 2|cj| ^ 



(A.IO) 



Substituting in Eq. ([5]) for the global entanglement we 
have 



AW- 



(A.ll) 



3 ij^j 



Separating the atom-mode terms and the mode-mode 
terms, and taking into account the symmetry of the two- 
qubit concurrence 



(A.12) 



it takes the form 



N 



where the basis states are \0j) and Using this we Q{W) — tt — 7 ( c^iPax) + E ^^(P 



first calculate and then its trace 



Am J 



\<\<ii<N 



trp2 = 1 + 2|c,f (I 



|co| 



l) , (A. 9) which leads to the result of Eqs. ® and (H 
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